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ABSTRACT 

We observed a sample of 20 representative Herbig Ae/Be stars and five A-type debris discs with PACS onboard of Herschel. The 
observations were done in spectroscopic mode, and cover far-IR lines of [O i], [C n], CO, CH + , H 2 and OH. We have a [O i] 63 fan 
detection rate of 100% for the Herbig Ae/Be and 0% for the debris discs. [Oi] 145 yum is only detected in 25 %, CO .7=18-17 in 45 
% (and less for higher J transitions) of the Herbig Ae/Be stars and for [C n] 157 yum, we often found spatially variable background 
contamination. We show the first detection of water in a Herbig Ae disc, HD 163296, which has a settled disc. Hydroxyl is detected 
as well in this disc. CH + , first seen in HD 100546, is now detected for the second time in a Herbig Ae star, HD 97048. We report 
fluxes for each line and use the observations as line diagnostics of the gas properties. Furthermore, we look for correlations between 
the strength of the emission lines and stellar or disc parameters, such as stellar luminosity, UV and X-ray flux, accretion rate, PAH 
band strength, and flaring. We find that the stellar UV flux is the dominant excitation mechanism of [Oi] 63, with the highest line 
fluxes found in those objects with a large amount of flaring and greatest PAH strength. Neither the amount of accretion nor the X-ray 
luminosity has an influence on the line strength. We find correlations between the line flux of [Oi] 63 and [Oi] 145, CO J = 18-17 
and [Oi] 6300 A, and between the continuum flux at 63 yum and at 1.3 mm, while we find weak correlations between the line flux of 
[O i] 63 and the PAH luminosity, the line flux of CO J = 3-2, the continuum flux at 63 yum, the stellar effective temperature and the 
Bry luminosity. Finally, we use a combination of the [Oi] 63 and 12 CO 7 = 2-1 line fluxes to obtain order of magnitude estimates of 
the disc gas masses, in agreement of the values we found from detailed modelling of 2 HAEBEs, HD 163296 and HD 169142. 

Key words. Stars - infrared, Astrochemistry, Line: identification, Protoplanetary discs 



1. Introduction 

Circumstellar discs around young stars are the sites of planet 
formation (e.g. Pollack et al. I1996I Alibert et al. 12005b . During 
the first 10 Myr, the initially gas-rich disc will evolve into first 
a transitional and then a debris disc, while dispersing its gas 
content. The understanding of this dispersal process and what 
favours/hinders it is a crucial part of the planet formation puz- 
zle, as the amount of gas present in a disc is crucial to determine 
whether gas giant planets can still be formed. Furthermore, the 
disc mass controls the migration of planetary bodies of all sizes, 
from gas giants to meter-sized planetesimals. Three components 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



need to be characterised well: the disc geometry, the dust, and 
the gas content. 

The disc geometry of young intermediate-mass stars, the 
Herbig Ae/Be stars (HAEBEs), is constrained through multi- 
wavelength imaging, interferometry, and radiative transfer mod- 
elling (e.g. Benisty et al. 120101) . Meeus et al. (2001) empiri- 
cally divided the HAEBE discs into group I (flared) and group 
II (flat). A general consensus exists that discs become flatter as 
dust grains grow and settle towards the midplane (Dullemond & 
Dominik 2004 ). Lately, several of the group I sources have been 
found to have an inner opacity hole in the disc (e.g. Grady et al. 
120071 120091 ), possibly due to a lack of small dust grains in the 
inner disc. In HD 100546, the gap may be caused by a planet 
(e.g. Bouwman et al. 120031 Tatulli et al. l20TT1 ). 
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In a study of 53 HAEBEs, 85% show a silicate emission fea- 
ture at 10 /mi with a variety in grain size and crystallinity, prov- 
ing the presence of warm small grains in these discs (Juhasz et al. 
120101 ). Polycyclic Aromatic Hydrocarbons (PAH) features were 
clearly detected in 70% of the sample, with a clear preference 
towards flared discs (Acke et al.|2010). PAHs located in the disc 
atmosphere are transiently excited by UV photons and are an im- 
portant heating source for the gas in the disc surface through the 
photo-electric effect. 

The study of gas properties is difficult as, in general, emis- 
sion lines are rather weak. Different gas species and transitions 
probe different regions in the disc: lines in the near- and mid-IR 
generally trace the inner disc (< 10-20 AU), while lines in the 
far-IR and mm mainly trace the outer disc. We refer to Carmona 
et al. (120101 1 for a discussion of different gas tracers, their loca- 
tion in the disc and observational characteristics. To understand 
the disc radial and vertical structure, it is necessary to observe 
several transitions of different species, as they arise under dif- 
ferent conditions (density, temperature, radiation field). H2 and 
CO lines are most often used, since they are the most abun- 
dant species present, with the canonical H2 to CO number ratio 
in the ISM being 10 4 . However, the detection of H2 in the IR 
has proven to be difficult because of its weak rotational and ro- 
vibrational transitions - it has only been detected in 3 HAEBEs. 
In a survey of 15 Herbig Ae/Be stars with CRIRES, Carmona 
et al. (1201 11 1 detected ro-vibrational transitions of H2 at 2.1218 
,1/m in only t wo ob jects: HD 97048 and H P 100546. Earlier, 
Bitner et al. ( 120081 ). Carmona et al. (2008), Martin-Zai'di et al. 
(20091 1201 Oi l searched for mid-IR pure rotational lines of H2 at 
17.035 jt/m in a sample of in total 20 HAEBEs; only two detec- 
tions were made, in AB Aur and HD 97048. In sharp contrast, 
the detection of CO, although much less abundant, is easier as its 
rotational/ro-vibrational lines are much stronger. CO is routinely 
detected in HA EBEs (e.g. Thi et al. 120011 Dent et al. I20051 . 
Lorenzetti et al. (2002) showed ISO/LWS observations of atomic 
and molecular lines in the far-IR for a sample of HAEBEs. They 
detected the fine-structure lines of [0 1] 63 & 145.5 fim and [C 11] 
at 157.7 fim. 

Despite the wealth of observations, it is still not clear 
how HAEBE discs dissipate with time. In the less massive T 
Tauri stars, disc dispersal is thought to be initiated by photo- 
evaporation, mainly due to ionising EUV (hv > 13.6 eV) pho- 
tons that first create a gap in the inner disc, which is subse- 
quently rapidly viscously accreted. In a next step, the outer disc 
is efficiently removed through a photo-evaporative disc wind 
(Alexander et al. 2006). However, Gorti et al. (2009) showed 
that UV can rapidly disperse the outer disc, where the bulk of 
the disc mass is located, thus setting the disc lifetime. Also X- 
rays are thought to play an important role in those discs (e.g. 
Ercolano et al. 120081 Owen et al., 120121 . And finally, also the 
accretion of a planet with a mass of a few Jupiter can play an 
important role in the dissipation of the disc. Which mechanism 
is ultimately dominating the dispersion process is not yet deter- 
mined. 

We present ESA Herschel Space Observatory (Pilbratt et al. 
2010 ) spectroscopy of 20 HAEBEs and 5 A-type debris discs, 
covering several transitions of abundant atoms and molecules 
that can be used as crucial tests of our understanding of disc 
physics and chemistry in the upper layers of the disc. The ob- 
servations cover a significant part of the disc surface that was 
not accessible before. Our observations are part of the Herschel 
Open Time Key Programme (OTKP) "GAS in Protoplanetary 
Systems" (GASPS; P.I. Dent; see Dent et al. l20T2l . With this 
paper, we want to obtain a better understanding of HAEBE discs 



by relating several gas tracers and excitation mechanisms with 
stellar and disc properties. What gas species are present in a 
HAEBE disc and at what temperatures? What is the physical and 
chemical structure of the disc chemistry? What is the dominant 
excitation mechanism of gas in HAEBE discs? 

In Sect.|2j we describe the sample and our methods to de- 
rive the stellar and disc parameters. In Sect. [3] we present the 
spectroscopy and the line detections. We discuss gas lines as a 
diagnostic tool in Sect. [4] and look for correlations between the 
observed line fluxes. We relate our detections and upper limits to 
stellar and disc parameters in Sect.|5] Finally, we round off with 
conclusions in Sect. [6] 

2. Targets 

The sample consists of 20 Herbig Ae/Be stars with spectral types 
ranging between B9.5 and F4, to which we will refer as Herbig 
Ae (HAe) stars. We do not include the more massive Herbig 
Be stars which are, in general, younger and have smaller discs 
and often an additional remnant envelope (e.g. Natta et al. 2000 
Verhoeff et al. 120121 . We also include HD 141569 A, an object 
that completely lacks a near-IR excess for A < 4.5 fim, attributed 
to inner disc clearing, but still has a substantial amount of pri- 
mordial gas; in this paper, we will call this a transitional disc. 
We are aware that several of our sources are also called (pre- 
)transitional discs in the literature, such as HD 100546 (Grady 
et al. 120031 ) and HD 135344B (Andrews et al. l20TTl ). which are 
observed to have a gap in their disc, but for the purpose of this 
paper, we include them in the HAe sample, since they still have 
a substantial near-IR excess and their total IR excess is much 
larger than that of HD 141569A. Besides the Herbig Ae sample, 
we include five debris discs around A-type stars with ages be- 
tween ~10 and 1200 Myr for comparison, as the Herbig Ae stars 
are seen as their precursors. We list the main stellar parameters 
in Table[T] The sample is representative for the known Herbig Ae 
stars: there are nine objects in Meeus group I and ten in group II 
(to which we will refer to as flaring and flat discs). Furthermore, 
we have a good coverage of T e g, age, stellar luminosity L„ and 
accretion rate. 49 Cet is the only debris disc in our sample for 
which gas was detected through CO observations with the JCMT 
(Dent et al. 120051) . Hughes et al. (120081) later resolved the CO gas 
emission using the SMA. 

For 3 stars in our sample, members of the GASPS team 
have performed detailed modeling of their discs with the ra- 
diative transfer code MCFOST (Pinte et al. [20061 120091 and 
the thermo-chemical code ProDiMo (Woitke et al. |2009al ): HD 
169142 (Meeus et al. [20 101), H D 100546 (Thi et al. ML§b and 
HD 163296 (Tilling et al. l20l2l . 

2.1. Stellar and disc properties 

To characterize the sample in a consistent way, we first aimed 
at determining the stellar component of the spectral energy dis- 
tribution (SED). In the next step, several parameters that can be 
important in the context of gas excitation in the disc were com- 
puted, namely, ultraviolet luminosities, infrared excesses, and 
accretion luminosities. 

We compiled a set of literature and catalogue stellar pa- 
rameters (effective temperatures, gravities and metallicities), 
and critically selected what we considered the best. Multi- 
wavelength photometry from different sources and ultravio- 
let spectra obtained by the International Ultraviolet Explore^ 

1 http://sdc.cab.inta-csic.es/ines/ 
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Table 1. Main stellar parameters of the sample. 



Star 


Alternative 


Sp. Type 


TeB 


logs. 


[M/H] 


Refs. 


d 




A v 


Age 




Name 




(K) 








(PC) 






(Myr) 


AB Aur 


HD 31293 


AO Ve 


9280 


4.00 


0.00 


1 


139.3+19.0 


33.0 ±9.2 


0.25 


5.0 ± 1.0 


HD 31648 


MWC 480 


A3-5 Ve 


8250 


4.00 


0.00 


2 


137.0 ± 26.2 


13.7 ±5.5 


0.16 


8.5 ±2.0 


HD 35187 




A2 Ve 


8915 


4.00 


0.00 


3 


114.2 ±32.4 


17.4 ± 10.6 


0.78 


9.0 ±2.0 


HD 36112 


MWC 758 


A5 IVe 


7750 


3.50 


-0.08° 


4 


279.3 ± 75.0 


33.7 ± 19.3 


0.16 


3.7 ±2.0 












-0.14* 


4 








3.5 ±2.0 


CQTau 


HD 36910 


F3 Ve 


6900 


4.35 


0.00 


5 


113.0 ±24.0 


3.4 ± 1.5 


1.40 


4.0 ±2.0 


HD 97048 


CU Cha 


AO Ve 


10000 


4.00 


-0.05" 


4 


158.5 ± 15.7 


30.7 ±6.1 


1.15 


6.5 ± 1.0 












-0.75 c 


4 








4.0 ±0.5 


HD 100453 




A9 Ve 


7400 


4.20 


+0.30 


6 


121.5 ± 9.7 


8.8 ± 1.4 


0.00 


> 10 












-0.02" 


4 








10.0 ±2.0 












-0.09* 


4 








9.0 ±2.0 


HD 100546 




B9 Ve 


10470 


3.50 


-0.08" 


3,4,7 


96.9 ± 4.0 


22.7 ± 1.9 


0.09 


- 












-13tf 


4 








3.8 ±0.5 


HD 104237 


DX Cha 


A4-5 Ve 


8550 


3.90 


+0.16 


8 


114.7 ±4.7 


28.8 ±2.4 


0.16 


5.5 + 0.5 


HD 135344 B 


SAO 206462 


F3-4 Ve 


6810 


4.40 


+0.14" 


10,4 


142.0 ± 27.0 


8.1 ± 3.1 


0.37 


10.0 ±2.0 












-0.08* 


4 








9.5 + 2.0 


HD 139614 




A7 Ve 


7400 


4.00 


-0.50 


6 


140 ± 42 


7.6 ±4.6 


0.00 


9.2 ±2.0 












-0.19° 


4 








11.5 ±2.0 












-0.27* 


4 








10.5 ± 2.0 


HD 141569 A 




B9.5 Ve 


10000 


4.28 


-0.50 


11 


116.1 ± 8.1 


29.6 ±4.2 


0.37 


4.7 ± 0.3 


HD 142527 




F6 Hie 


6250 


3.62 


0.00 


12 


233.1 ±56.2 


33.2 ± 16.9 


0.59 


2.0 ±0.5 


HD 142666 




A8 Ve 


7500 


4.30 


+0.20 


6 


145 ± 43 


13.5 ±8.0 


0.93 


9.0 ±2.0 


HD 144668 


HR 5999 


A7IVe 


7925 


4.00 


0.00 


13 


162.9 ± 15.3 


50.8 ± 9.5 


0.47 


2.8 ± 1.0 


HD 150193 


MWC 863 


A2 IVe 


8970 


3.99 


0.00 


2 


216.5 ±76.0 


48.7 ± 38.0 


1.55 


3.8 ±2.0 


KK Oph A 




A6 Ve 


8000 


4.00 


0.00 


16 


260 


13.7 


1.60 


8.0 ±2.0 


KK Oph B 




G5 Ve 


5750 


4.50 


0.00 


17 


260 


2.1 


2.80 


8.0 ±2.0 


51 Oph 


HD 158643 


B9.5 Hie 


10250 


3.57 


+0.10 


2 


124.4 ± 3.7 


285.0 ± 17.0 


0.37 


0.7 + 0.5 


HD 163296 


MWC 275 


Al Ve 


9250 


4.07 


+0.20 


14 


1 18.6 ± 11.1 


33.1 ±6.2 


0.47 


5.5 ±0.5 


HD 169142 


MWC 925 


A7-8 Ve 


7500 


4.00 


-0.50 


15 


145 ± 43 


9.4 ±5.6 


0.00 


7.7 ±2.0 


49 Cet 


HD 9672 


A4 V 


9500 


4.30 


+0.10 


2 


59.4 ± 1.0 


21.0 ±0.7 


0.22 


o Q+6.1 A1+1I9 

°- y -2A' UI -46 


HD 32297 




AO V 


9520 


4.15 


0.00 


18 


112.4 ± 10.8 


10.9 ±2.1 


0.62 




HR 1998 


HD 38678, f Lep 


A2IV-V 


8500 


4.27 


-0.76 


9, 19 


21.6 ±0.1 


14.0 ±0.1 


0.00 


1250 ± 250 


HR 4796 A 


HD 109573 A 


AO Ve 


9750 


4.32 


0.00 


9 


72.8 ± 1.8 


23.4 ± 1.1 


0.00 


10.0 ±2.0 


HD 158352 


HR 6507 


A7 V 


7500 


3.85 


0.00 


9 


59.6 ±0.9 


17.7 ±0.6 


0.00 


1000 + 200 



Notes: Quantities in italics are assigned, (a) These metallicities are weighted averages of all the elemental abundances listed in Tables 2, 3 and 4 
of Acke & Waelkens (2004), see Appendix A of Montesinos et al. (2009) for details; (b/c/d) Metallicities assumed to be in the same proportion as 
the species [Fe i]/([Si n])/[Fe n], respectively (see Table 3 of Acke & Waelkens, 2004). Refs.: (1) Woitke et al. (GASPS) (in prep), (2) Montesinos 
et al. (2009), (3) Manoj et al. (2006), (4) Acke & Waelkens (2004), (5) Testi et al. (2003), (6) Guimaraes et al (2006), (7) Levenhagen et al. (2006), 
(8) Fumel & Bohm (201 1), (9) Allende-Prieto et al. (1999), (10) Miiller et al. (201 1), (1 1) Merin et al. (2004), (12) Verhoeff et al. (201 1), (13) van 
Boekel et al. (2005), (14) Tilling et al. (GASPS) (2012), (15) Meeus et al. (GASPS) (2010), (16) Herbig (2005), (17) Carmona et al. (2007), (18) 
Torres et al. (2006), (19) Gray (2006). Distances are from the revised parallaxes by van Leeuwen (2007) except for HD 135344 B (Miiller et al. 
201 1) and HD 139614, HD 142666 and HD 169142 (van Boekel et al. 2005). 



were used to construct the spectral energy distribution of the 
star-plus-disc systems. A specific model photosphere for each 
star was extracted, or computed by interpolation, from the 
grid of PHOENIX/GAIA models (Brott & Hauschildt (20051). A 
Rayleigh- Jeans extension up to mm wavelengths was applied to 
the original models. The model photosphere was reddened with 
several values of E(B - V) (Ry = 3.1) and normalized to the flux 
at V (0.55 yum), until the best fit to the optical photometry was 
achieved. In some cases, the photospheric model fits the ultra- 
violet spectrum fairly well, implying only stellar photospheric 
contribution to that spectral range, whereas in other cases an ex- 
cess, attributed to accretion processes, was apparent. 

The photospheric, infrared, and ultraviolet fluxes (required 
to estimate the stellar luminosity, L IR /L„ and L uv /L ), were all 
computed from the dereddened data or model. The photospheric 
flux was estimated by integrating the model photosphere; the 
infrared flux was computed by fitting a spline to the observed 
data set, integrating from the wavelength Ao where the fit started 
separating from the photospheric model up to 200 fim; the UV 



flux between 1 150-2430 A was used to compute Luv- Accretion 
rates were estimated from the observed photometric excesses 
in the Balmer discontinuity. These were modelled assuming a 
magnetospheric accretion geometry, following the procedure by 
Mendigutfa et al. (1201 1 bb . Upper limits on the accretion rates 
are provided for most stars, given that the Balmer excesses are 
negligible. 

In Table[T]we give some of the basic parameters of the stars, 
namely, identifications, spectral types, effective temperatures, 
gravities, metallicities and the corresponding references (cols. 
1-7), distances (col. 8), and two parameters derived from those 
and the SEDs, namely, stellar luminosities and ages (cols. 9 and 
10). In Table|2]we list the disc groups according to Meeus et al. 
(2001) or the evolutionary status of the stars, fractional infrared 
luminosities, ultraviolet luminosities and accretion luminosities 
(cols. 2-5). Stars with two values of L\jy/L in col. 4 show two 
different emissivity levels in IUE spectra obtained in different 
epochs; the numbers corresponding to the low and high states 
are given. 
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Table 2. Derived properties of the sample. Disc group classifi- 
cation from Meeus et al. ( 1200 11 1 and Acke et al. (2010). For the 
derivation of L\ R /L t , L V y/L q and log L. dcc /L Q , see the descrip- 
tion in Sect l2.ll 



Star 


Disc 
group 


Lm/U 




log L , cc /La 


AB Aur 


I 


0.76 


4.63 


< 0.55 


HD 31648 


11 


0.46 


0.75 


< 0.19 


HD 35187 


11 


0.12 


2.23 


1.13 


HD 36112 


1 


0.66 


1.32 


< -0.81 


CQTau 


II 


1.01 


0.19-0.3 


< -1.16 


HD 97048 


I 


0.39 


7.69 


< 0.95 


HD 100453 


I 


0.62 


0.29 


< -0.91 


HD 100546 


1 


0.56 


7.22 


1.62 


HD 104237 


II 


0.32 


1.54 


0.87 


HD 135344 B 


I 


0.56 


0.11 a 


-0.22 


HD 139614 


I 


0.39 


0.39 


< -1.12 


HD 141569 A 


II/TO 


0.009 


6.83 


< 0.70 


HD 142527 


I 


0.98 


0.15 11 


< -1.04 


HD 142666 


11 


0.33 


0.37-0.68 


0.85 


HD 144668 


11 


0.51 


1.55-2.94 


< 0.23 


HD 150193 


II 


0.48 


8.53 


1.32 


KK Oph A+B 


II 


2.01 


2.35 


1.70 


51 Oph 


11 


0.03 


71.32 


< 0.37 


HD 163296 


11 


0.29 


3.21-5.58 


-0.04 


HD 169142 


I 


0.42 


0.45 


< 0.05 


49 Cet 


Debris 


6.0 xlO" 4 


2.96 


< -0.11 


HD 32297 


Debris 


0.003 


1.91 a 




HR 1998 


Debris 


6.6 x 10" 5 


1.25 




HR 4796 A 


Debris 


0.003 


5.59 




HD 158352 


Debris 


1.6 X 10" 4 


0.64 





Notes: Disc groups I and II refer to Meeus et al. (2001) classification, 
"TO" stands for "Transition Object", (a) UV fluxes measured on the 
model photosphere (no UV observations available). 



The age estimations were done by placing the values of 
log T e ff - log L t /L Q for each star on an HR diagram containing 
tracks and isochrones computed for its particular metallicity. For 
some stars this parameter is unknown, therefore a solar abun- 
dance was assumed, which introduces an uncertainty which is 
difficult to estimate. Metallicity is an important parameter to take 
into account when determining ages using this procedure since 
the position of a set of tracks and isochrones in the HR diagram 
changes substantially with metallicity. The evolutionary tracks 
and isochrones for a scaled solar mixture from the Yonsei-Yale 
group (Yi et al. 2001 ) - Y 2 in their notation - have been used in 
this work. From the Y 2 set, the isochrones with Z = 0.02 (solar) 
have been used for those stars with measured metallicities [M/H] 
between -0.10 and +0.10 and for the stars for which no abun- 
dance determinations are available. The remaining metallicities 
are treated with the isochrones in the grid whose value of Z is 
closer to 0.02 x 10 [M/H] . 

There is one star (HD 32297) for which a determination of 
the age was impossible; its position falls below the main se- 
quence in a Z = 0.02 HR diagram, therefore, it is quite likely 
that this star has subsolar abundance. The evolutionary stage of 
49 Cet is interesting and was studied in detail in Montesinos et 
al. (|2009). Two ages, corresponding to PMS and MS isochrones, 
are listed in Table Q] 

2.1.1. KK Oph 

Special attention had to be paid to KK Oph. This object is a close 
binary separated by 1.61 arcsec (300 AU at our adopted distance, 



Table 3. Ranges and lines targeted with the PACS spectrometer. 



Set 


Observed Range 
A (/mi) 


Species 


Transition 


Wavelength 

(jum) 


A 


62.68-63.68 


[Oi] 


s Pi ^ 'Pi 


63.184 






o-H 2 


8l8 — > 7o7 


63.324 




188.77- 190.30 


DCO+ 


/ = 22-21 


189.570 


B 


71.90-73.05 


o-HoO 


7(>7 — * 6 1 6 


71.946 






CH + 


J = 5-4 


72.14 






CO 


J = 36-35 


72.843 




144.0- 146.1 


p-H 2 


4i3 -> 3 22 


144.518 






CO 


J = 18-17 


144.784 






[Oi] 


3 P -> 3 Pi 


145.525 


C 


78.55 - 79.45 


o-H 2 


4 23 -> 3 a 


78.741 






OH 


1/2 - - 3/2 


79.11 






OH 


1/2 +- 3/2 


79.18 






CO 


J = 33-32 


79.360 




157.1 - 158.9 


[Cn] 


2 Pl/2 ~* ~Pl/2 


157.741 






p-H 2 


hi -> 4 04 


158.309 


D 


89.45 - 90.50 


p-H 2 


3 22 — > 2u 


89.988 






CH+ 


J = 4-3 


90.02 






CO 


J = 29-28 


90.163 




178.9-181.0 


o-H 2 


2 12 — > loi 


179.527 






CH+ 


J = 2-1 


179.610 






o-H 2 


2 2 i — > 2u 


180.488 



see discussion below; Leinert et al., [T997). the hot component 
(A) is a Herbig Ae star and the cool component (B) is a T Tauri 
star, both of them actively accreting (Herbig, 120051 and refer- 
ences therein). The SED shows an infrared excess from ~ 1 //m 
onwards. All the available photometry corresponds to both com- 
ponents, therefore it is not straightforward to estimate parameters 
for each star from these data alone. In addition, the system is 
highly variable (Hillenbrand et al. |1992t Herbst & Shevchenko, 
1999; Eiroa et al. 120011 Oudmaijer et al. 120011 adding further 
complications to the analysis. 

Herbig (2005 ) determined an spectral type A6 V for KK Oph 
A, and Carmona et al. (2007) - who only studied the secondary 
component - found a spectral type G5 V for KK Oph B. In 
both cases the authors made use of high resolution spectroscopy 
which allowed them to separate the spectra of each component. 
Effective temperatures of 8000 and 5750 K have been assigned 
to the stars according to their spectral types. To estimate lumi- 
nosities from the available photometry we have followed the ap- 
proach by Carmona et al. d20071 l where each star is affected by 
a different extinction, adopting Ay — 1.6 and 2.8 mag for com- 
ponents A and B, respectively. We have worked on the set of 
optical and near-IR photometry provided by Hillenbrand et al. 
(|1992| l, and assumed that the contribution to the flux at U and 
B from the cool component was negligible. The distance to this 
system is uncertain, with published values of 170-200 pc. We 
find these distances too low to obtain values for the luminosities 
of the stars that are in agreement with their spectral type and lu- 
minosity class; a lower limit of ~ 260 pc provides sensible results 
and has been adopted in this work, which is the best we could 
derive with the available data set. 

3. Herschel PACS spectroscopy 

We obtained PACS (Poglitsch et al. 2010) spectroscopy 
in both line and range modes (PacsLineSpec, 1669 s and 
PacsRangeSpec, 5150 s). The observation identifiers (obsids) 
can be found in Table [A] 1 of Appendix A. In a later stage, we 
obtained deeper range scans to confirm tentative detections by 
doubling the integration time for 7 sources. All the observations 



4 



G. Meeus et al.: Herschel GASPS observations of HAEBEs 



were carried out in ChopNod mode, in order to remove the emis- 
sion of the telescope and background. PACS is an IFU with 25 
spaxels, 9. "4 on each side. Due to the characteristics of the PSF 
at 60yum, only ~70 % of a point source flux falls in the cen- 
tral spaxel, with a decrease towards longer wavelengths, down 
to 45% at 180 jum. 

The spectroscopic data were reduced with the official re- 
lease version 8.0.1 of the Herschel Interactive Processing 
Environment (HIPE; Ott 2010 ), using standard tasks provided in 
HIPE. These include bad pixel flagging; chop on/off subtraction; 
spectral response function division; rebinning with oversample 
= 2 and upsample =1, corresponding to the native resolution of 
the instrument; spectral flatfielding and finally averaging of the 
two nod positions. In order to conserve the best signal and not to 
introduce additional noise, we only extracted the central spaxel, 
and corrected for the flux loss with an aperture correction pro- 
vided by the PACS instrument team ('pointSourceLossCorrec- 
tion.py'). This is only possible when the source is well-centered. 
If that was not the case (e.g. for HD 142666), several spaxels 
needed to be taken into consideration to obtain the full target's 
emission. 

The main lines targeted are the fine structure lines of [Oi], 
[C n] and the molecular lines of CO, OH, CH + and H 2 0. In to- 
tal we observed eight spectral regions; the details of the tran- 
sitions can be found in Table[3] The spectral resolution varies 
between 3400 (shortest wavelengths) and 1100 (longest wave- 
lengths), which are equivalent to ~ 88 kms~' at 60 pm, and 177 
km s _1 at 190 pm. In our sample, we do not resolve the emission 
lines as our objects do not have such high-velocity components. 
The line flux sensitivity is of the order a few 10~ 18 W/m 2 . The 
lines were identified by manual inspection of automated gaus- 
sian fits at the expected line position, taking the instrumental 
resolution as the expected FWHM of the line. We extracted the 
fluxes of the detected lines using a gaussian fit to the emission 
lines with a first-order polynomial to the continuum. We used 
the RMS on the continuum (excluding the line) to derive a lcr 
error on the line by integrating a gaussian with height equal to 
the continuum RMS and width of the instrumental FWHM. This 
approach is necessary as HIPE currently does not deliver errors 
for the spectra. In case of a non-detection, we give a 3cr upper 
limit, also calculated from the continuum RMS. The absolute 
flux calibration error given by the PACS instrument team is cur- 
rently <15 %. The measured line fluxes - in case of detection - 
or their upper limits are listed in Tables |4] and [5] Several atomic 
and molecular lines were observed, as described in the following 
paragraphs. 



3.1. Oxygen fine structure lines 



The third most abundant element in the interstellar medium is 
oxygen. Its fine structure line at 63.2 pm is by far the strongest 
line observed in our spectra (see Table [4]). It is detected in all 
Herbig Ae/Be stars, and absent in the debris discs. For five of our 
objects, the line flux is larger than 200 x 10~ 18 W/m 2 , while our 
faintest detection is 20 x 10~ 18 W/m 2 . The other fine structure 
line [O i] 145 pm is also one of the strongest lines in our spectra, 
however it is only detected in five objects (25% of the HAEBEs). 
In Figs. [T] and [2] we show the spectra centered on [Oi] 63 and 
[O i] 145 pm for the whole sample. 



Table 4. Atomic line strengths in units of 10~ 18 W/m 2 with lcr 
continuum RMS between brackets in case of a detection or 3cr 
upper limits. For HD 158352, not all settings were observed, 
therefore we write 'n.a.' (not available) when the spectrum is 
lacking. ": an absorption feature is observed in the spectrum, that 
is caused by the subtraction of stronger emission in the chop-off 
position. 



Star 




roii 


rc Hi 


A (yjm) 


63.18 


145.53 
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^ 1 o 
< /.O 


^ Ci 1 

< y. / 


I 1 U jj 1 o / 


O I A Q1 

jZ.O (4.oJ 


^ A f. 


< J.y 
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Al O 


< 4.0 


^ fx A 
< 0.4 


HD 139614 


44.5 (6.1) 


< 4.7 


<8.3 


HD 141569A 


245.3 (4.8) 


24.9 (1.4) 


11.4(2.1) 


HD 142527 


52.3 (3.8) 


< 11.6 


<28.7 


HD 142666 


18.9 (3.1) 


<4.7 


<9.0 


HD 144668 


140.2 (4.3) 


<5.8 


<5.4 


HD 150193 


24.6 (2.8) 


< 6.6 


<5.3 


KK Oph 


172.8 (5.0) 


6.2(1.3) 


9.2(1.2) 


51 Oph 


53.3 (2.5) 


<5.1 


<6.8 


HD 163296 


208.4 (4.2) 


<4.0 


ABS" 


HD 169142 


91.5 (4.4) 


<3.6 


< 8.0 


49 Cet 


< 10.0 


< 6.3 


<9.1 


HD 32297 


<7.4 


<4.2 


<7.1 


HR 1998 


<7.0 


<5.1 


<6.6 


HR 4796A 


< 6.3 


<3.8 


<5.4 


HD 158352 


< 8.2 


n.a. 


<3.8 



3.2. Carbon fine structure line 

When detected, the [Cn] 157.7 //m line can be strong - more 
than 100 x 10~ 18 W/m 2 . However, it is only seen in six ob- 
jects (30% of the HAEBEs; see Fig. [3]l - these are the same 
objects for which [Oi] 145 fim was also detected (see Table HJ, 
apart from HD 36112. The low detection rate for [Cn] is sur- 
prising, given the high detection rate (83%) in ISO/LWS spec- 
tra reported by Lorenzetti al. (2002). This can be attributed to 
a difference in aperture: 80 arcsec for LWS versus 9 arcsec for 
PACS. Moreover, we noticed that the line can also be present 
in the off-source chop positions (in a spatially variable amount), 
contaminating our spectra. In two cases, the dominating emis- 
sion is present in all on and off-source spaxels, so that our chop- 
off subtracted spectra even show the feature in absorption as the 
chop-off position shows stronger emission (HD 142527 and HD 
163296, see Fig . [3]> . The interpretation of the [C n] emission line 
is complex. Besides originating in the disc, it could also form in 
the remaining envelope, or simply in cloud material in the line 
of sight. A detailed analysis of background [C n] emission is be- 
yond the scope of this paper, but a dedicated study is underway 
(Pantin et al., in preparation). 



3.3. Carbon monoxide 

In our PACS ranges, we cover four transitions of the CO 
molecule: 7=36-35 (72.85 ^m), 7=33-32 (79.36 pm), 7=29-28 
(90.16/mi) and 7=18-17 (144.78 pm). These are all mid to 
higher 7 transitions, with E llp between ~ 950 and 3700 K. The 
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Fig. 1. The [O i] 63 fim lines for the entire sample. The line is seen in emission in all the HAEBEs, while absent in the more evolved 
debris disc objects 49 Cet, HD 32297, HR 1998, HR 4796A and HD 158352 (bottom row). 



highest 7 transitions in our settings (7=36-35 and 7=33-32) are 
not detected in any star of our sample. In Fig. |4] we show the 
region around 90 /mi, covering the CO 7=29-28 transition. This 
CO line is only clearly seen in AB Aur and HD 100546, with 
a tentative detection for HD 97048. For the lowest 7 observable 
(18-17, at 144.78 /on) we see many more detections in our spec- 
tra (see Pig . [2j : it is detected in nine objects, and there is one 
tentative detection, for HD 169142. The strongest CO lines are 
observed in AB Aur, HD 97048 and HD 100546. HD 141569A is 
the only star for which we have a clear (more than 5<x) detection 
of [Oi] 145 fim, but no detection of CO at 144.78 //m, showing 



that both species trace different excitation conditions and chem- 
istry (atomic versus molecular). We will discuss this further in 
Sect.1431 



3.4. Hydroxyl 

With our settings, we only cover the OH doublet at 
79.1 1/79. 18 /mi. In Fig. [5] we can see evidence for the doublet 
in HD 100546 and HD 163296; however, the only 3 cr detection 
(of one line of this doublet) is found in AB Aur and HD 163296 
(see Table [6}. Sturm et al. (2010) detected several OH lines at 
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Table 6. Molecular line strengths of H2O, hydroxyl and CH + for the sources with a least one detection in these lines. Line flux in 
units of 1(T 18 W/m 2 , lcr continuum RMS between brackets in case of a detection or, in case of a non-detection, 3<x upper limits. 
With ':' we indicate that the feature is between 2 and 3<x. At 90 and 179.5 fim, there is a possible blend of CH + and H2O. For the 
water lines, we also give the upper energy level in K. 



Star 


o-H 2 


o-H 2 


CH + 


0-H2O 


OH 


OH 


CH7 p-H 2 


CH + / o-H 2 


0-H2O 


A (fim) 


63.32 
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72.14 
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79.11 


79.18 
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Fig. 2. The settings around 145 /urn. The lines of CO J=18-17 at 
144.8 fim and [Oi] 145 fim are only clearly detected in a three 
objects. 



far-IR wavelengths (53-200 fim) in the SED-mode PACS spec- 
tra of HD 100546. 



3.5. Water 

The only star in our sample with convincing evidence for water 
is HD 163296. The measured line fluxes are listed in Table|6] 
While the feature at 63.32 fim is seen in HD 163296 with 3<x 
confidence, we list it as a tentative detection given the spurious 
absorption feature next to it (in HD 31648, a potential feature 
of water at 63.32 fim is redshifted by 0.009 fim so we also con- 
sider it a tentative detection). Also in our deeper range scans at 
71.946 and 78.74 fim we see evidence for water in HD 163296 
(see Fig. [5}. The features are close to 3<x, however, the fact that 
we do see emission lines at the positions where water lines are 
predicted to be present further strengthens the detection of water 
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Fig. 3. The [Cn] line at 157.7 fim. The line is only clearly de- 
tected in three objects, more weakly in three others and is ap- 
pears in absorption in HD 142527 and HD 163296, due to sub- 
traction of the chop-off spectrum that has stronger emission. 



in HD 163296. The analysis of the water i n HP 163296 will be 
presented in a separate paper (Meeus et al. I2Q121 >. 

In two other objects, HD 97048 and HD 100546, we do see 
several emission lines at the position of water. In Fig.|4]we show 
the region around 90 micron, where the transition of para-H20 
can appear. However, when we zoom in on these stars in the re- 
gions where other water lines are expected to be present, we get 
a different picture; in Fig.|5]we show the spectra covering water 
lines. At 90.00 fim (see Fig. [5} we do see an emission peak for 
both stars, which could be the paraTTiO line (at 89.988 fim), but 
it is a blend with CH + at 90.02 fim. At other positions of water, 
which are not blended with CH + , we do not see any emission 
line: 71.948, 78.741 and 78.93, 180.488/im). 

To summarise, the only times we detect an emission feature 
at the position of a water line in these two objects, is in a blend 
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Fig. 5. Comparison of HD 100546 (top spectrum), HD 97048 (middle) and HD 163296 (bottom) at 72.6 and 79.0, 90.0 and 180.0 
/im. We indicate the position of the lines of CO, CH + , OH and H2O. 



with CH + . We can conclude that there are no H2O emission lines 
detected in these two objects at the sensitivity of our observa- 
tions. Therefore, we do not confirm the detection of water in HD 
100546 reported by Sturm et al. (2010) as also seen in a later 
improved reduction of the original data (Sturm, priv. comm.). 

3.6. CH + 



Thi et al d201 II ) reported the first detection of CH + in a Herbig 
Ae/Be disc, HD 100546. We detect the feature of CH + in both 
HD 97048 and HD100546 at 90.00 ^m. For HD100546 we have 
detections of two more lines, at 72.14 and 179.61 jjm. The fea- 
tures at 90.0 and 179.6/im are also at the position of water fea- 
tures, but given the lack of other water features, they can be at- 
tributed to CH + . 



4. Gas lines as tracers of the conditions in the disc 

The far-IR lines observed with PACS form in different regions 
of the disc. The exact locations vary with geometry, i.e. flaring, 
inner holes, gaps, but we will discuss here for simplicity the gen- 
eral case of a continuous flaring disc around a Herbig Ae star. 

The [C n] line depends strongly on the irradiation of the star, 
especially UV photons shortward of 1200 A (Pinte et al. 120101 
Kamp et al. 201 1). The emission originates foremost in the upper 
tenuous layers of the disc (low critical density) where UV pho- 
tons can penetrate. The [0 1] 63 and [0 1] 145 //m form deeper in 



the disc where the atomic oxygen abundance is still high. Most 
of their emission comes from 10-100 AU as the temperatures be- 
yond 100 AU are generally too low to excite those lines (Kamp 
et al. 1201 Oi l. The high excitation water lines (i^per > 400 K) 
form mostly in the surface layers of the hot water reservoir in- 
side the snow line (15 AU for an effective temperature of 10500 
K, moving inwards for cooler stars while keeping the disc struc- 
ture constant). The low excitation lines (^upper < 200 K) form 
in a thin layer beyond the snow line where water can be pho- 
todesorbed from the icy grains into the gas phase (Cernicharo et 
al. 120091 Woitke et al. l2009bb . The exception is the 89.988 jum 
water line with an upper level energy of ~ 300 K which forms 
across the snow line (for TW Hya, a T Tauri disc - Kamp et al. in 
preparation). Bruderer et al. (2012) modelled the CO ladder in 
HD 100546. They found that the high J lines of CO can only be 
reproduced by a warm atmosphere in which T gas is much higher 
than Td U st- The low J lines of CO (observed in the mm) trace the 
outer disc (at several 100 AU radial distance), while the mid to 
high J lines observed in the far-IR originate at distance of several 
tens of AU. The highest J lines of CO form mostly in the very 
inner disk (typically inside a few AU, or at the rim of transition 
discs). 

The fundamental ro-vibrational CO band (Av = 1) at 4.7 /zm 
band, tracing the terrestrial planet-forming region is routinely 
observed in HAEBEs (e.g. Brittain et al. [2007). The bands are 
rotationally excited up to high J(> 30), with r rot between 900 
and 2500 K (van der Plas et al. l20"l2l . If the gas is not in LTE, 
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Table 5. CO line strengths in units of 10" 18 W/m 2 , l<x contin- 
uum RMS between brackets in case of a detection or, in case of 
a non-detection, 3<x upper limits. With ':' we indicate that the 
feature is between 2 and 3 <x. For HD 158352, not all settings 
were observed, therefore we write 'n.a.' (not available) when the 
spectrum is lacking. 



H ? CH + COJ29-28 



then the vibrational temperature, TVa,, can depart from T lot , when 
UV fluorescence causes super-thermal level populations. This is 
observed in several UV bright HAEBEs where 7\,;b > 5000 K: 
HD 97048 and HD 1 00546 have T vib > 6000 K, wh ile T ml ~ 
1000 K (Brittain et al. 120071 van der Plas et al.EJn). Also, in 
group I discs r rot < r v ;b, while in group II discs, TVa, ^ r iot (van 
der Plas etal. 

Furthermore, the line profile suggests CO depletion in the 
innermost regions of HAEBE discs (van der Plas et al. 2009), 
with group I clearing a larger radius (r,„ ~ 10 AU) than group 
II discs (r,„ ~ 1 AU; van der Plas et al. 120121 . The transitional 
disc HD 14 1569 A stands out for having a low T ml ~ 250 K, 
while its TVa, (~ 5600 K) is in a similar range of the hottest CO 
observed in HAEBEs, attributed to UV fluorescence (Brittain et 
al. l2007l l. Furthermore, Goto et al. d20061 l showed that this disc 
has an inner clearing in CO up to a radius of 1 1 AU, comparable 
to the group I discs. 

12 CO lines observed in the millimeter come from low / tran- 
sitions of optically thick CO located in the outer disc surface. 
These pure rotational transitions of cold CO (Av = 0) are rou- 
tin ely de tected in HAEBE discs (e.g. Pietu et al. 120031 Dent et 
al. [2005 ). Earlier, the existence of Keplerian rotation in discs 
was confirmed with mm interferometry of CO lines (Koerner et 
al. 1 1993b . Furthermore, as the lines are optically thick, a simple 
model of the line profile allows for an estimate of the outer disc 
radius and even inclination (e.g. Dent et al. 120051 Panic et al. 
2008). 



Star 


CO J36-35 


CO J33-32 


CO J29-28 


CO J18-17 


A (yum) 


72.85 


79.36 


90.16 


144.78 


AB Aur 


< 37.9 


< 33.4 


31.0 (4.3) 


26.6 (8.8) 


HD 31648 


<9.6 


<16.3 


<9.6 


6.9 (2.3) 


HD 35187 


< 21.0 


<28.9 


< 8.4 


< 5.5 


HD 36112 


< 30.2 


< 17.6 


< 17.7 


6.6(1.1) 


CQTau 


< 15.4 


< 25.4 


<9.9 


5.9(1.5) 


HD 97048 


< 21.3 


< 16.6 


15.3: (6.0) 


28.5 (2.6) 


HD 100453 


< 13.9 


< 19.5 


< 9.6 


<8.3 


HD 100546 


< 28.7 


< 38.8 


49.3 (16.1) 


58.0(3.3) 


HD 104237 


< 19.6 


< 16.1 


<9.8 


<5.6 


HD 135344B 


< 13.0 


< 21.7 


< 8.5 


< 8.4 


HD 139614 


< 16.4 


< 19.9 


<9.0 


<6.9 


HD 141569A 


< 12.3 


< 10.6 


< 6.6 


<3.4 


HD 142527 


< 30.9 


< 23.7 


< 12.4 


<9.8 


HD 142666 


< 16.3 


< 18.8 


< 12.6 


< 6.2 


HD 144668 


< 11.9 


< 12.1 


< 8.2 


5.6(1.3) 


HD 150193 


< 14.8 


< 13.1 


<7.2 


< 6.2 


KK Oph 


<6.5 


< 9.4 


< 7.7 


6.9(1.3) 


51 Oph 


<9.0 


< 12.4 


<6.6 


<4.0 


HD 163296 


<8.0 


<10.7 


<7.9 


11.6(1.5) 


HD 169142 


< 10.0 


< 10.9 


< 10.0 


5.1: (1.7) 


49 Cet 


< 13.6 


< 18.6 


<9.8 


<6.7 


HD 32297 


< 11.4 


< 16.9 


<9.7 


<4.3 


HR 1998 


< 12.6 


< 13.5 


< 6.6 


<4.8 


HR 4796A 


< 12.8 


< 17.1 


<7.6 


<4.1 


HD 158352 


n.a. 


<9.2 


n.a. 


n.a. 
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Fig. 4. The spectra around 90 micron, covering P-H2O at 89.988, 
CH + at 90.02 and CO at 90. 163 pm. The only sources with 
clearly detected lines are AB Aur (CO), HD 97048 and HD 
100546 (CH + and CO). HD 100546 is scaled by a factor 1/4, 
so that the emission lines would still fit in the plot window. 



Woitke et al. (2010) calculated a grid of disc models with 
the thermo-chemical radiation code ProDiMo. This model grid, 
called "Disc Evolution with Neat Theory" is a useful tool to 
derive statistically meaningful dependencies on stellar and disc 
properties. Kamp et al. (1201 li ) used the model grid to derive line 
diagnostics that are relevant for the PACS observations. We will 
refer to these diagnostics in our discussion below. Our sample 
includes several objects that stood out in earlier papers, in terms 
of detections of H2, CO and/or OH, which can be attributed to a 
high level of UV fluorescence; we will relate these results also 
to our new observations. We will now discuss the results pre- 
sented in the previous sections in the context of our current un- 
derstanding of these discs. The following sections present our 
interpretation of observational correlations and their meaning in 
the context of disc structure and evolution. 

In order to remove the bias that can be introduced by the 
distance of the stars, we scaled our data to a distance of 140 
pc (to ease comparison with objects in Taurus and predictions 
of the model grid). The scaled data are all the PACS line and 
continuum fluxes, the I2 CO J = 3-2 and 2-1 line fluxes, and the 
mm continuum fluxes. 

The relations between parameters are analysed with their 
corresponding 'p-values' (see Table [7]), which gives the prob- 
ability that the two variables considered are not correlated. Two 
parameters will be classified as 'correlated' if one or more of the 
three (Spearman, Kendall and Cox-Hazard) differently obtained 
p-values are not larger than 1%, and as 'tentatively correlated' 
when 1< p < 5% (e.g. Bross 1197 ll ). When p > 5 %, the pa- 
rameters are classified as 'not correlated', as their /^-values are 
similar to ^-values derived for randomly generated samples. The 
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Table 7. Probability p (in percentage) that the two parameters (x,y) under consideration are not correlated, calculated with several 
statistical methods: Spearman's, Kendall's and Cox-Hazard's. Under 'Result' we indicate what we can derive from these statistical 
tests. 



Parameters 

X 


y 


p Spearman 


p Kendall 


p Cox-Hazard 


Result 


Linear fit: y 
a 


=axx+b 
b 


log [0 1] 63yum 


log 


[Oi] 145yum 


0.07 


0.06 


<0.1 


Correlated 


1.2 + 0.21 


1.8 ±3.3 


log [0 1] 63yum 


log 


CO 718-17 


0.06 


0.03 


<0.1 


Correlated 


0.64 ± 0.07 


-7.1 ± 1.1 


log Luv 


log [0 1] 63/jm 


0.33 


0.32 


- 


Correlated 


0.72 ± 0.20 


-16.1 ± 0.1 


log T eS 


log 


Luv 


0.05 


<0.1 


<0.1 


Correlated 


9.3 ± 0.2 


-36.5 ± 4.4 


log L[oi]6300A 


log 


[O i] 63yimi 


0.88 


0.45 


- 


Correlated 


0.83 ±0.21 


-12.4 ± 1=0.9 


log F 63/Im 


log 


F 1.3mm 


0.19 


0.09 


0.01 


Correlated 


0.72 ±0.1 


-0.20 ± 0.2 


log F 63jUm 


log 


[0 1] 63yimi 


4.6 


4.4 


2.7 


Correlated? 






log CO J 3-2 


log 


[0 1] 63yimi 


3.7 


2.4 


- 


Correlated? 






log r eff 


log 


[O i] 63yum 


4.7 


3.5 


1.6 


Correlated? 






log Lp A H6.2/im 


log 


[O i] 63yum 


6.6 


4.4 




Correlated? 






log L Brr 


log 


[0 1] 63yimi 


4.6 


4.4 




Correlated? 






log CO J 2-1 


log 


[0 1] 63yum 


54 


31 




Not correlated 






log L, 


log 


[O i] 63yum 


7.2 


9.2 


18 


Not correlated 






logLx 


log 


[0 1] 63yimi 


84 


71 




Not correlated 






log F 13mm 


log 


[O i] 63yimi 


34 


27 


44 


Not correlated 






log L Acc (BaDis) 


log [O i] 63yum 


84 


83 




Not correlated 






log Fir/F, 


log [O i] 63yum 


21 


27 


47 


Not correlated 






slope b 


log 


[O i] 63yimi 


43 


38 


72 


Not correlated 







p-values and linear fits provided in Table [7] take into account 
that several of our datasets include upper/lower limits instead 
of detections, as they were derived with the ASURV package 
(Feige lson & Nelson [1985] Isobe et al. 1 19861 and Lavalley et al. 
1992 ), that was specifically designed to deal with censored data. 
We used the Spearman's partial correlation technique to quantify 
the influence of the common distance parameter on the proba- 
bility of false correlation and found that this is negligible for our 
sample - the p-values considering the distances, or random val- 
ues instead, are practically equal. The absence of influence of the 
distances on the correlations most probably comes from the rela- 
tively narrow range covered by this parameter in our sample. We 
excluded 51 Oph in the correlation test with Luv and [Oi] 63, as 
it is an outlier when comparing its extremely high Luv to the rest 
of the sample (see Tableland Fig. [8j, and an enigmatic object 
(e.g. van den Ancker et al. 2Q0Tb. 

4.1. Oxygen fine structure lines 

In Fig.[6| we plot the strength of the [O i] 63 /mi line as a function 
of the continuum flux at 63 pm. The variables are weakly corre- 
lated (see Table|7|. The three sources with the highest line fluxes 
also have the highest continuum fluxes: AB Aur, HD 97048 and 
HD 100546, to which we will refer to as 'the bright three'. These 
are also the only HAEBEs in which H2 emission has been de- 
tected in the IR (see Sect.[T]). 

The ratios of the fine structure line fluxes of [0 1] 63, 
[Oi] 145 and [Cn] 157 pm are diagnostics of the excitation 
mechanism (e.g. Kaufman et al. 1999). Unfortunately, for most 
of the sources we only obtained upper limits for one or more 
of these lines. We show the line flux of [Oi] 63 as a function of 
[Oi] 145 in Fig.|6| and find a clear correlation (see Table[7]). We 
find line ratios of [Oi]63/[Oi] 145 between 10 and 30. These 
ratios are not compatible with predictions of the PDR model in 
Tielens & Hollenbach (1985 ) for optically thick lines with T gas < 
200 K. Our line ratios (median of 24) are in agreement with pre- 
dictions from the model grid, which gives a median line ratio of 
25 (Kamp et al. 1201 11 1. These authors conclude, based on those 
disc models, that the line ratio is not sensitive to the average 



oxygen gas temperature (for 50 < 7 gas < 500 K), but instead 
correlates with the gas to dust ratio. 

4.2. Ionised carbon fine structure line 

The [Cn] line flux is very sensitive to the UV radiation field, 
and the line is mostly optically thin. Unfortunately, most of our 
sources are background contaminated, as [C 11] is also detected 
in off-source positions, in a variable amount, depending on the 
location. For the few sources with solid, non-[Cn] background 
contaminated detections, we find line ratios of [0 1] 63/[C 11] 157 
between 10 and 30. 

4.3. Carbon monoxide 

Freeze-out of CO on grains is not expected in the disc of an 
A-type star (e.g. Panic et al. 2009), so that the strength of the 
low-/ 12 CO lines can be used to obtain a lower limit on the 
cold gas mass. The disc size can be derived from the 12 CO flux 
and profile. In Fig. [7] we plot the [Oi] 63 line flux as a function 
of the I2 CO 7=2-1 and 7=3-2 line strengths (data from Dent 
et al [20051 Panic et al. 120091 Isella et al. I20T01 Oberg et al. 
l20101l20TTl and our own data, see Appendix |Bj). We do not find 
a clear correlation with the 7=2-1 transition, but we do find a 
weak correlation with the 7=3-2 transition (see Table[7]). 

Kamp et al. ( |20TTT ) showed that the ratio of [Oi] 63/ 12 CO 
7=2-1 can be used to derive the gas mass in the disc to an or- 
der of magnitude. The idea is that this ratio is determined by the 
average gas temperature in the disc. If the [Oi]63 line is opti- 
cally thin, the line flux will mainly depend on the gas mass and 
average [Oi] temperature (Woitke et al. 120 101 ). Once the temper- 
ature is known, the line flux of [0 1] 63 can thus be related to the 
disc gas mass. In our sample, we have nine sources for which 
the 12 CO 7=2-1 line flux is known. We calculated the ratio for 
those sources, and found that log ([Oi] 63/ 12 CO 7=2-1) falls be- 
tween 2.5 and 3.5. This means that we have a similar average 
gas temperature in all cases. We apply the relation for this ratio 
range between log [0 1] 63 and the gas mass, derived by Kamp et 
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Table 8. Line fluxes of 12 CO 7=2-1, log of the line ratios and 
derived gas masses. 



Object 


12 CO 7=2-1 
(1CT' 8 W/m 2 ) 


log([Oi] 63/ 
CO 7=2-1) 


M gils 
(Mo) 


HD 31648 


0.169 


2.75 


6.5 xl0" J 


HD 35187 


0.026 


3.10 


2.5 xl0~ 3 


HD 36112 


0.099 


2.58 


0.24 xl0~ 3 


CQTau 


0.024 


3.30 


4.4xl0- 3 


HD 135344 B 


0.080 


2.78 


2.2 xlO- 3 


HD 139614 


0.054 


2.92 


2.2 xl0~ 3 


HD 142527 


0.160 


2.52 


0.66xl0- 3 


HD 142666 


0.052 


2.57 


0.62xl0- 3 


HD 163296 


0.053 


3.59 


25.4x10-3 


HD 169142 


0.164 


2.75 


5.3 xlO" 3 


49 Cet 


0.015 


2.84 


< 3.3 xlO" 3 



al. d201 11 1 to obtain an estimate of the disc mass. The results are 
shown in Table [£] We derive M gas between 0.24 and 25 x 10" 3 
M Q . These values are of course only indicative; for a more ac- 
curate estimate, a full model of all the available observations 
needs to be done for each disc. The masses are consistent with 
the estimates derived from a detailed modelling of HD 163296 
(M gas ~ 15 - 120 x IO" 3 M Q ; Tilling et al. |20T2b and HD 169142 
(M ~ 3 - 7 x IO" 3 M ; Meeus et al. l2010b . 



4.4. Hydroxy I 

Although less abundant than H2 and CO, hydroxyl (OH) is also 
an important molecule, as it plays a central role in the forma- 
tion/destruction of H 2 0, H 2 and [0 1]. Mandell et al. ( 12008b were 
the first to detect ro-vibrational transitions of warm OH (at 3.0- 
3.7 fim) in two HAEBEs, AB Aur and HD 361 12. They derived 
a rotational temperature of 650-800 K, and argue that fluores- 
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cent excitation is responsible for the emission of OH located 
in the disc surface layer. Fedele et al. (201 lj) also searched for 
OH in 11 HAEBEs with CRIRES, and detected it in 4 sources 
with spectral types between B5 and Al; none of those objects 
are in our sample. They find that objects with an OH detection 
tend to be Meeus group I sources. Recently, several transitions of 
OH around 3 /jm were detected in HD 100546 (Liskowsky et al. 
I20T21 . In our PACS spectra, we only detect OH in HD 163296 
(and a tentative detection in AB Aur) - but we only cover one 
doublet, which is not the strongest in the far-IR, so it might be 
we are not sensitive enough. Indeed, Sturm et al. (2010) detected 
several OH lines in their full SED range mode PACS spectra of 
HD 100546, from which the 84 fim doublet is the strongest. In 
our spectrum at 79 fim, we do see a indications for the OH dou- 
blet, but it is not a 3<x detection (see Fig.[5]l. 

4.5. Water 

Water was not yet reported to be detected in a Herbig Ae/Be 
disc despite several searches in the near- and mid-IR (e.g. 
Pontoppidan et al. 2010| Fedele et al. 1201 II) . However, in the 
study by Pontoppidan et al. d20T0i l. HD163296 shows a H 2 
emission line at 29.85 /mi, but in that paper, water was only con- 
firmed to be present when detected both at 15.17 and 17.22 jjm, 
with at least 3.5cr confidence. 

Thi & Bik d20031 ) showed that the ratio H 2 0/OH declines 
when the ratio of the intensity of the UV field over the density in- 
creases. Thus in lower density regions with a lot of UV radiation, 
the amount of water expected is low, compared to OH. Indeed, 
Fedele et al. J2.01 1 b conclude that, if water vapour is present, it 
must be located in deeper, colder layers of the disc than where 
OH is found; the disc atmosphere is depleted in water molecules. 

We detect at least one water line and have evidence for sev- 
eral others in HD 163296, a group II source. Tilling et al. (12012b 
modelled the disc of HD 163296 based on our earlier, shallower, 
range scans, and showed that the disc is mostly settled, which 
results in slightly warmer dust and increased line flux. This fact, 
together with the rather high UV luminosity, can probably ex- 
plain the water detections in this disc. 

Kamp et al. ( 1201 11 1 found that strong dust settling will in- 
crease the water abundance in the disc surface. The reason is 
complex (we refer the reader to Sect. 5.3 of Kamp et al.), but 
the main idea is that there is an efficient cold water formation 
route in these discs. It is interesting to note that HD 100546, the 
source that has the highest UV flux and is richest in other, strong 
lines ([Oi]63, CO and CH + ) does not show evidence of warm 
water. HD 100546 is a group I source, thought to have its in- 
ner disc cleared (e.g. Bouwman et al. 12003 1 Benisty et al. 2010). 
Woitke et al. d2009bl l showed that water lines originate in 3 dis- 
tinct regions: 1) a deep midplane behind the inner disc wall, up 
to 10 AU, hosting most of the water vapour; 2) a midplane re- 
gion between 20 and 150 AU where water freezes out and there 
is a small amount of cold water vapour; and 3) a warm water 
layer between 1 and 50 AU higher up in the disc. In HD 100546, 
region 1 and part of region 3 are missing, so that the amount of 
gas phase water is much smaller. 

Furthermore, while in the inner disc of HD 163296 the den- 
sity is too low for water to form (from OH + H 2 ) to balance the 
fast photodissociation (see e.g. Thi & Bik 2005), water can sur- 
vive in the inner 10 AU of the warm atmosphere. In contrast, in 
the UV strong star HD 100546, even at 30 AU the UV field is 
too strong for water to survive. This, in combination with the 
greater amount of settling in the HD 163296 disc, might ex- 



plain the absence of detectable H 2 emission in the disc of HD 
100546 while it is detected in HD 163296. 



4.6. CH + 

The formation of CH + is controlled by the gas-phase reaction 
C + + H 2 — > CH + + H, which has an activation energy of 4500 
K. Therefore, CH + not only traces the presence of H 2 but also 
the presence of hot gas. In our sample, we have IR detections 
of molecular hydrogen in 3 targets: AB Aur, HD 97048 and HD 
100546. In two of these objects, we also detect CH + , suggest- 
ing that their formation and excitation mechanisms are indeed 
related. For a more in-depth discussion of CH + in HD 100546 
we refer to Thi et al. d20TT1) . 



5. Correlations with stellar and disc parameters 

In the next paragraphs, we will look for correlations between the 
observed [Oi]63 line fluxes and the properties of the objects. 
For this purpose, we did not include the debris discs, as their 
discs have a very different nature. Besides, in the debris discs, 
we did not detect any [0 1] 63 line, so that we would only be able 
to compare upper limits. 



5.1. The influence of r eff , UV and X-ray luminosity 

We searched for a correlation between stellar parameters and the 
[0 1] 63 line flux. We did not find a trend with age nor with stellar 
luminosity. In Fig. [8] we show the relation between the line flux 
of [Oi]63 and the effective temperature of the stars. Both ap- 
pear uncorrelated, until T e g reaches 10000 K, when the [Oi] 63 
flux increases dramatically for a few sources. Our /^-values (see 
Table |7]i are inconclusive, therefore, the significance of a possi- 
ble correlation cannot be established from our statistical analy- 
sis. 

The UV and X-ray photons play an important role in the 
chemistry and temperature balance of protoplanetary discs. For 
HAEBEs, the UV photons are important in heating the disc, 
through their absorption by PAHs and the subsequent photoelec- 
tric effect. In Fig. [8] we show the relation between the line flux 
of [0 1] 63 yum, the UV luminosity calculated from IUE spec- 
tra (see Sect.[2i, and the X-ray luminosity (data mainly found 
in Hubrig et al. 120091 see Table C.l for a full list). There is 



a clear correlation between the [0 1] 63 flux and the UV lumi- 
nosity (see Table |7J, as reported earlier for a limited sample in 
Pinte et al. (2010). This could be related to an increase in OH 
photo-dissocation in the disc surface and/or to a more efficient 
photoelectric heating of the gas by PAHs in those sources with a 
higher UV luminosity. 

X-ray photons can ionise atoms and molecules. The X-ray 
fluxes observed in HAEBEs (log Lx = 28-30) are on average 
lower when compared to those of the lower-mass T Tauri stars 
(TTS; log L x = 29-32). Aresu et al. J20TTT) found from theoreti- 
cal modeling of discs a correlation between the X-ray luminosity 
and the [0 1] 63 line flux for X-ray luminosities above 10 30 erg/s. 
Below that value, the gas temperature in the region where the 
[0 1] 63 line forms is dominated by UV heating while above that 
value, X-rays provide an additional heating source, thereby in- 
creasing the total line flux. Since all objects in our sample are 
below this Lx threshold, it is not surprising that we do not see 
a correlation with the X-ray luminosity. Furthermore, X-rays in 
HAEBEs are softer than in TTS, so that they cannot penetrate as 
deep in the discs as in T Tauri discs. A dedicated study will use 
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spectral X-ray properties to interpret the observed PACS spectra 
(Giidel et al., private communication). 

5.2. Relation with accretion rate 

In Fig. [9] we plot the line flux of [O i] 63 jum as a function of 
L acc derived from the excess in the Balmer discontinuity. There 
is no trend visible, just more scattering at higher L acc . However, 
we should point to the difficulties that lie in an accurate determi- 
nation of the accretion rate, that will introduce additional scatter 



in the values. We also plot the [Oi] 63 line flux against the Bry 
luminosity (data from Garcfa-Lopez et al.|2006 and Donehew et 
al. I201 II) . Here we see a tentative correlation with L^ry Given 
that we do not see a clear correlation between the accretion rate 
and the [O i] 63 line flux, we can conclude that the accretion 
is not an important contributor to the excitation of [O i] 63 in 
HAEBEs. This confirms the findings of an earlier study of a few 
HAEBEs by Pinte et al. d20T0b . that the emission from HAEBE 
discs can be explained by photospheric heating alone. Indeed, in 
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Fig. [9] bottom, we plot Luv vs. T e g, and see that both are very 
well correlated. This means that the bulk of the UV luminosity 
in HAEBEs is photospheric, rather than originating in accretion, 
as is often observed in the cooler T Tauri stars (e.g. Yang et al. 
I20T21 . 

5.3. Relation with PAH and [Ol]6300 A 



PAHs are important for the heating of the disc through the photo- 
electric effect. In HAEBEs, the PAH luminosity, Lpah/L* is ob- 
served to reach up to 9 xlO 3 . In Fig. 10 we show the relation 
between the PAH luminosity and the [0 1] 63 line fluxes (data 
from Acke et al.|2004bjand Keller et al. EoUBl The [Oi] 63 flux 
and the PAH flux weakly correlate with each other (see Table|7]i. 

HD 141569A is the only star for which we saw a firm de- 
tection of [0 1] 145 pm, but no detection of CO at 144.8 pm. The 
[0 1] 63 over [0 1] 145 ratio is smaller than 10, while it is around 
20 in ABAur, HD 97048 and HD 100546. Moreover, Brittain 
et al. d2007l > showed from CO 4.7 pm modelling that T mt is of 
the order of 250 K, while T ml in HD 97048 and HD 100546 is 
much higher, > 1000 K. This difference cannot be attributed to a 
lower UV luminosity as it is rather similar (6.83 in HD 141569A 
vs. 7.69 in HD 97048 and 7.22 L G in HD 100546). However, a 
lower PAH luminosity is observed in HD 141569A (a factor 10 
less; Acke et al. 1201 Oi l, so there is less heat contribution to the 



disc - there might also be an intrinsic difference in PAH abun- 
dance. 

The excitation mechanism of [0 1] 6300 A is still not well 
determined, it could be either thermal or non-thermal (e.g. flu- 
orescence or OH photodissociation). The emission line is com- 
monly found in HAEBEs (e.g. Corcoran & Ray 1998). Acke et 
al. ([2005 ) derived from a simple model of the spectroscopically 
resolved [Oi] 6300 A line that the emission originates from the 
disc surface, assuming that it is non-thermally excited by the 
UV photo-dissociation of OH molecules in the disc surface. It 
is more often found in flared than in self-shadowed discs (Acke 
et al. 12005b and traces the disc between 0.1 and 50 AU (van der 
Plas et al. 120091 . In the bottom plot of Fig. [10] w e plot [0 1] 63 
against [Oi] 6300 A (data from Acke et al. 120051 van der Plas 



et al. 2008 and Mendigutfa et al. 2011a). Here we see a corre- 
lation that could be explained by the excitation mechanism for 
both lines being related to UV photons. We already showed that 
there is no obvious correlation between the accretion rate and the 
[0 1] 63 line strength. The special properties of the 'bright three', 
which have a high UV luminosity as witnessed by the high line 
flux of [0 1] 63 and PAH luminosity, thus must mainly be be- 
cause of their high effective temperature, when compared to the 
rest of the sample. The only other objects in our sample with a 
temperature around 10000 K, are HD 14 1569 A (a transitional 
disc) and 51 Oph, a special case with a sharp drop in the SED at 
longer (> 20 micron) wavelengths, which can be attributed to a 
compact disc in which the dust has settled towards the midplane. 

5.4. Relation with disc properties 

The continuum flux at 1.3 mm is often used to derive a mini- 
mum dust mass of the disc, under the assumption that the dust 
emission is optically thin (Beckwith et al. 1 1990b . In Fig. 1 1 
we plot the [0 1] 63 flux as a function of the continuum flux 
at 1.3 mm (our mm data collected in the literature are listed in 
Appendix C.2 we also use our own unpublished SMA data, see 
Appendix BJ. We did not find that these variables are correlated 
(see Table[7Ji. On the other hand, the continuum flux at 63 /mi 



and at 1.3 mm are strongly correlated (see Fig. Hi. 

In Fig.[l2j we plot the line flux of [0 1] 63 pm as a function 
of the total IR excess, a proxy for the amount of dust continuum 
observed. We do not find a correlation. This is likely because the 
IR continuum is rather a tracer of the dust disc, and the scale 
height of the gas disc can be higher than that of the dust, as 
is already observed in a few HAEBEs (e.g. van der Plas et al. 
2009). HD 14 1569 A has a transitional disc, with a much lower 
Fi R /F„ than the rest of the sample (see Table [2]), indicating that 
the disc is barely flaring and perhaps has already dissipated much 
of its disc material. All these diagnostics lead to confirm that the 
disc of HD 141569A is much different (inner disc mostly cleared 
out from dust) from that of HD 97048 and HD 100546, which 
still have flaring gas-rich discs. Modelling of the HD 14 1569 A 
disc with ProDiMo will help us to better understand these dif- 
ferences in terms of excitation mechanisms, abundance of gas 
and disc structure (Thi, in preparation). Our data suggest that, 
in HD 14 1569 A, the CO is located deeper in the disc (closer to 
the dust) where it can be thermalised and/or shielded from pho- 
todissociation by direct UV photons. 

The slope b of the far-IR to mm SED, where F^ ~ A h , can be 
related to the size of the dust grains radiating at mm wavelengths. 
However, grain size is not the only factor influencing the slope, 
also the composition (e.g. amount of carbon) and grain shape can 
be important factors. Acke et al. (2004a) showed that the SED 
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Fig. 12. Top: [Oi]63 versus amount of IR excess. We do not 
find a correlation between these variables. The 2 objects to the 
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disc, respectively. Bottom: [O i] 63 versus the slope b of the SED 
at far-IR to mm wavelengths. The [O i] 63 flux also does not cor- 
relate with the SED slope. Diamonds: group I sources, asterix: 
group II sources. 



far-IR to mm slope is related to the SED group: self-shadowed 
discs (group II) have on average shallower slopes than their flar- 
ing counterparts (group I). In our sample, we do not see a corre- 
lation between the [O i] 63 line strength and the SED slope (see 
Fig. ITT). 



5.5. Non-detections in debris discs 

In our sample we have five debris discs, for which the [O i] 63 
line was not detected. We obtain 3cr upper limits for the 
[Oi]63 line flux (at their respective distances) ~ 6-10 xl0~ 18 
W/m 2 . This contrasts with the young debris disc f3 Pic, where 
[O i] 63 fim and [C n] 158 fim emission lines were detected with 
Herschel/PACS (Brandeker et al. 2012). These authors give an 
[O i] 63 urn line flux = 13.2 xlO -18 W/m 2 for p Pic (at a distance 
of 19 pc), what would not have been detected at the distance 
of our debris discs. The only exception is HR 1998, at a dis- 
tance of 22 pc, for which we have an upper limit ~ 5 xl0~ 18 
W/m 2 when scaled to 19 pc, almost a factor 3 lower than the 
P Pic detection, p Pic appears to be a special debris disc that 
is relatively rich in gas, originating from the ongoing vaporisa- 
tion of dust through grain-grain collisions, comet evaporation, 



and/or photodesorption of grain surfaces (Lagrange et al. 1998; 
Czechowski & Mann [20071 Chen et al. |2"00"7I 

6. Conclusions 

In this paper, we studied the gas content with Herschel PACS 
spectroscopy for a sample of 20 Herbig Ae stars and five A-type 
debris discs, that can be summarised as follows: 

1 . We detect the [O i] 63 yum line in all the Herbig Ae stars of 
our sample, while it is absent in the debris discs, confirming 
the lack of a large amount of gas in these discs. The [O i] 63 
line is by far the strongest line observed in our spectra, next 
in strength (if detected) are [O i] 145 and [C n]; they are only 
detected in 5 (25%) and 6 (30%) sources, respectively. 

2. The CO mid to high J transitions (18-17 and 29-28) are only 
detected in 9 (45%) and 2 (10%) objects, respectively. The 
highest J (33-32 and 36-35) CO lines covered in our spec- 
tra are not seen at all in our sample. The three detections of 
CO 7=29-28 are in the three strongest UV emitting objects, 
AB Aur, HD 97048 and HD 100546, revealing the need for 
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a large amount of UV photons for this line to become visi- 
ble. Interesting in this respect is the transitional disc of HD 
141569 A, where we did not detect CO 7=18-17, but did de- 
tect a strong line of [Oi] 145. This cannot be attributed to 
a difference in UV luminosity but rather to significant inner 
disc clearing, and to a more tenous disc. 

3. We detect two lines of CH + in HD 100546, and also detect 
CH + for the first time in HD 97048, only the second Herbig 
Ae star in which it is detected. 

4. Hydroxyl and H2O are important ingredients of the disc 
chemistry. However, we found water and OH in only one 
object, HD 163296, which has a settled disc. The previ- 
ous detection of H2O, announced by Sturm et al. (1201 Oi l in 
HD 100546 cannot be confirmed. The misidentification was 
caused by a blend with the CH + line, often present at the 
same wavelength as H2O. The non-detection of H2O in most 
sources is in agreement with findings of Pontoppidan et al. 
(120 10]) and Fedele et al. (120111) . who also did not detect water 
at IR wavelengths, despite dedicated surveys. 

We correlated the strength of the [0 1] 63 fim line with stellar 
parameters, as well as disc properties. We can summarise our 
findings as follows: 

1 . The [0 1] 63 line flux correlates weakly with the continuum 
flux at 63 /im. The line flux ratios of [Oi] 63/ [Oi] 145 and 
[0 1] 63/ [C n] are between 10 and 30. 

2. We found that three of our sources, AB Aur, HD 97048 and 
HD 100546, have very strong [0 1] 63 line fluxes, when com- 
pared to the rest of the sample. These three sources have 
group I discs and have the highest T e ff values in the sam- 
ple, and thus have more stellar UV flux. Indeed, we see a 
correlation between the total (stellar + accretion) UV lumi- 
nosity and the strength of the [Oi] 63 line. We do not see a 
correlation with the X-ray luminosity, which is rather low in 
our sample of HAEBE stars. 

3. We did not find a correlation between the accretion rate es- 
timated from the B aimer discontinuity, and a tentative one 
with the Bry line. This shows that accretion is not the main 
driver of the [0 1] 63 excitation in HAEBEs. The bulk of the 
UV luminosity is photospheric rather than from accretion. 

4. Sources with high [0 1] 63 fluxes also have high PAH lumi- 
nosity, which can both be related to their high UV fluxes. We 
also see a correlation with the luminosity of the [0 1] 6300 A 
line. 

5. The disc geometry (flat versus flared) does not uniquely 
determine the strength of the [0 1] 63 line flux. The three 
strongest lines are observed in flared discs, but once these 
sources are excluded, there is no significant difference in line 
strength observed between the group I and II discs. 

6. We found a strong correlation between the continuum at 
63//m and at 1.3 mm. There is no correlation between the 
[0 1] 63 line strength and the strength of the dust continuum 
at 1.3 mm. We also did not find a correlation with the slope 
of the far-IR to mm SED, nor with the IR excess. 

7. We see a weak correlation with the strength of 12 CO 7 = 3-2 
line. Based on the line ratio [0 1] 63/ 12 CO 7=2-1, we can de- 
rive an estimate of the gas mass present in the disc. We found 
M gas between 0.25 and 25 x 10~ 3 M , consistent with the 
estimates derived from a detailed modelling of HD 163296 
(M gas ~ 15 - 120 x lO -3 M ; Tilling et al. 120121 and HD 
169142 (M gas ~ 3 - 6.5 x 10~ 3 M ; Meeus et al. l20T0t . 

A picture emerges for the protoplanetary discs around 
Herbig Ae/Be stars where the stellar UV flux is the main param- 



eter controlling the strength of the [0 1] 63 line, which is formed 
just below the disc surface. An increased amount of settling can 
enhance the line flux for those species (such as water or OH) 
that are formed deeper in the disc, where the density is higher. 
We plan to follow-up on this study with detailed modelling of 
a few key objects: AB Aur & HD 97048 and HD 135344 B & 
HD 142527 (group I, high and low UV, respectively), HD 163296 
(group II), HD141569 A (transitional disc), and finally the enig- 
matic compact disc of 5 1 Oph. Our modelling results will further 
aid in the understanding of the chemistry and physical processes 
present in Herbig Ae/Be discs. 
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Appendix A: PACS observation identifications 

In Table |Aj we show the Herschel observation identification 
numbers (obsids) of our observations. Several stars were ob- 
served twice in the range mode, to obtain a deeper observation 
These are also indicated in the table. In Figs 
we show all the stars at 72, 79 and 180 /im. 



A.l A.2andA.3 



Appendix B: Additional millimeter observations 

Interferometric observations of HD 32297, HD 35187, HD 
142666, HD 144668, HD 158352, HR 1998, and KK Oph were 
carried out in five tracks with the Sub-Millimeter Array (SMA; 
Ho et al. 2004) in the compact configuration from November 
2010 to October 2011. The observations were carried out at a 
central frequency of 223.9 GHz, with upper and lower band- 
passes of 4 GHz bandwidth and the center of each bandpass 
offset from the central frequency by 2.5 GHz. The continuum 
was sampled at 3.25 MHz, and a 104 MHz regions was set aside 
for 0.203 MHz resolution observations of the CO J=2-l line 
(230.5379700 GHz) line. The compact array observations in- 
cluded seven or eight antennas, with baselines from 10-70 m. 
Each object was observed for ~2.5 hours total integration time 
in good conditions (zenith T225GHZ 0.06 - 0.25) with system tem- 
peratures of -100-230 K . 

For gain calibration, we interleaved 5 minute observations of 
close quasars between 20 minute observations of the targets. We 
combined observations of two or three targets per track. We used 
60-90 minute observations of bright quasars for bandpass cali- 
bration, and 20 minute observations of available planets for flux 
amplitude calibration. Observations of flux and bandpass cali- 
brators were carried out before or after our targets were avail- 
able. We flagged and calibrated the data using standard routines 
in the facility IDL package MIR. We carried out baseline based 
phase calibration, finding rms phase errors of ~ 10-20°. Based 



17 



G. Meeus et al.: Herschel GASPS observations of HAEBEs 



Table A.l. Overview of the obsids that were observed. (D) 
means that it was a deeper observation than our regular settings; 
(Dl) for settings at 79/158 //m, (D2) for settings at 72/145 /mi 
and 79/158 //m and (D3) for all 3 range settings. 



CH + COJ36-35 



Star 


Line Spec 


Range Spec 


HD 9672 


1342188424 


1342188423 


AB Aur 


1342191355 


1342191354 






1342226000 (D2) 


HD 31648 


1342226002 


1342226003 


HD 32297 


1342217849 


1342217850 


HD 35187 


1342226900 


1342226901 


HD 36112 


1342227635 


1342227636 


HD 36910 


1342227638 


1342227637 


HR 1998 


1342226192 


1342226193 


HD 97048 


1342188436 


1342188435 


HD 100453 


1342203059 


1342203058 




1342212228 


1342203442 


HD 100546 


1342188438 


1342188437 


HD 104237 


1342212234 


1342212233 


HR 4796A 


1342199242 


1342199243 






1342234488 (Dl) 


HD 135344B 


1342190370 


1342190369 


HD 139614 


1342191300 


1342191299 


HD 141569 


1342190376 


1342190375 






1342204340 (D3) 


HD 142527 


1342216173 


1342216172 


HD 142666 


1342214224 


1342214225 


HD 144668 


1342192146 


1342216200 


HD 150193 


1342216625 


1342216626 


KK Oph 


1342192148 


1342192149 






1342228205 (D2) 


HD 158352 


1342190377 


1342227800 (Dl) 


HD 158643 


1342217821 


1342217822 


HD 163296 


1342192161 


1342192160 






1342229742 (D3) 


HD 169142 


1342186310 


1342186309 






1342215676 (D2) 



Table B.l. 

SMA. 



i2r 



CO 2-1 and 12 CO 3-2 line fluxes observed with 



Star 


12 CO 2-1 Error 


12 CO 3-2 Error 




(Jy km/s) 


(Jy km/s) 


HD 35187 


< 0.55 




HD 142666 


<0.62 




HD 144668 


confused 


confused 


KK Oph 


< 1.08 




HR 1998 


< 1.17 




HD158352 


< 0.48 





on variations in the measured fluxes of our gain calibrators, the 
flux calibration has an uncertainty of ~15%. 

To measure the 1.3 mm continuum flux, we combined wide- 
band continuum channels from all observations and both side- 
bands for each object. We used the MIRIAD command uvfit to 
the fit a point source to the observations. In order to confirm that 
sources are unresolved in our observations, we carried out stan- 
dard Fourier inversion, CLEAN deconvolution and image recon- 
struction using natural weighting with the facility reduction tool 
MIRIAD. The typical synthesized beam of the combined obser- 
vations is ~ 3'.' 5 x 2'.' 8, with an rms of 0.3 - 1.5 mJy/beam in the 
continuum. We separately carried out the same imaging process 
on the higher spectral resolution line data, binning in 1.0 km/s 



as 



C 

<3 * 



AB Aur 



CQTau 



HD1 04237 



HD1 42527 



KK Oph 



49 Cet 



HD31648 



HD97048 




HD135344B 



HD142666 



51 Oph 



HD32297 



HD35187 



HD1 00453 



HD139614 



HD144668 



HD36112 



HD1 00546 



HD141569 



HD150193 



HD1 63296 



HR1998 



HD169142 



HR4796A 



72.0 72.5 73.0 72.0 72.5 73.0 72.0 72.5 73.0 72.0 72.5 73.0 

Wavelength (urn) 



Fig. A.l. The region around 72 micron. 
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Fig. A.2. The region around 79 micron. 



velocity channels centered on the CO J=2— >1 line. The beam- 
size is slightly smaller than in the continuum, with rms in each 
channel of 0.4 to 2 Jy/beam. We list the line fluxes in Table B.l 
and the 1.3 mm continuum fluxes in Table 
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O 
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Fig. A.3. The region around 180 micron. With a vertical line we 
indicate the position of CH + . 



Table C.l. X-ray luminosities collected from the literature, ref- 
erences listed: (l) Telles chi et al. d2007b . (2 ) Grad y et al. d2010l i, 
(3) Stelze r et al. d2004b . (4 ) Colli ns et al. ( 120091 (5 ) Stelz er et 
al. (120061 (6) Gr ady et al. (120091 , (7) S telzer e t al. d2009l) , (8) 
Berghoefer et al. (fT996t , (9) Grady et al. d2007l > and (10) Stelzer 
& Neuhauser J2000l. 



Star 




Instrument 


ref. 










AB Aur 


29.71 


XMM 


(i) 


HD 31648 


29.30 


Chandra 


(2) 


HD 97048 


29.58 


XMM 


(3) 


HD 100453 


28.82 


Chandra 


(4) 


HD 100546 


28.93 


Chandra 


(5) 


HD 104237 


30.11 


Chandra 


(5) 


HD 135344 B 


29.66 


ROSAT 


(6) 


HD 141569A 


< 28.1 


Chandra 


(5) 


HD 144668 


28.3 


Chandra 


(7) 


HD 150193 


29.64 


Chandra 


(5) 


51 Oph 


< 28.98 


ROSAT 


(8) 


HD 163296 


29.6 


Chandra 


(5) 


HD 169142 


29.1 


Chandra 


(9) 


HR 4796 A 


29.38 


ROSAT 


(10) 



Table C.l. Continuum fluxes at 1.3mm collected from the liter- 
ature; references listed: (1) Henning et al. ( 11994b . (2) Mannings 
et al. ( fl997l l. (3) Chapillon et al. (g fglli, (4) Meeus etal. (2 002 
12003]) ). (5) S ylveste r et al. ( fl996l (6) Acke et al. d2004al) . (7) 
Hughes et al. d2008l 



Table B.2. Continuum fluxes at 1.2 mm, observed with 
MAMBO and at 1.3 mm, observed with SMA. 



Star 


F 1.2mm 

(mJy) 


Error 


Fl .3mm 

(mJy) 


Error 


HD 32297 


3.14 


0.822 


3.10 


0.74 


HD 35187 


33.960 


1.001 


28.95 


0.85 


HD 141569 A 


4.785 


0.507 






HD 142666 


102.500 


2.022 


99.0 


4.3 


HD 144668 






34.3 


0.9 


HD 151093 


47.860 


2.091 






KK Oph 


36.160 


2.260 


24.5 


4.3 


51 Oph 


5.007 


0.599 






HR 4796 A 






<9.3 




HD 158352 






< 0.7 5 





Star Fi.3 mm Error Instrument ref. 
(mJy) (mJy) 



AB Aur 


103 


20 


IRAM 


(1) 


HD 31648 


360 


20 


OVRO 


(2) 


HD 36112 


56 


1 


IRAM 


(3) 


CQ Tau 


162 


2 


IRAM 


(3) 


HD 97048 


450 


30 


IRAM 


(1) 


HD 100453 


265 


21 


SIMBA/SEST 


(4) 


HD 100546 


470 


20 


IRAM 


(1) 


HD 104237 


90 


20 


IRAM 


(1) 


HD 135344B 


142 


19 


UKT14/JCMT 


(5) 


HD 139614 


242 


15 


UKT14/JCMT 


(5) 


HD 142527 


1190 


30 


ATCA 


(6) 


HD 163296 


780 


30 


OVRO 


(2) 


HD 169142 


197 


15 


UKT14/JCMT 


(5) 


49 Cet 


2.3 


0.6 


SMA 


(7) 



In addition, HD 32297, HD 35187, HD 142666, KK Oph, 
HD 141569, HD 150193, and HD 158643 were observed for 
continuum emission at 1 .2 mm using the MAMB02 bolometer 
array (Kreysa et al. 1998) on the IRAM 30m telescope at Pico 
Veleta, Spain. Observations were conducted during the Nov. 
2008 bolometer pool. Zenith opacity for our observations was 
typically ~0.2 - 0.3, and observations were carried out to a tar- 
get lcr sensitivity of 1 mJy, typically 20 minutes on source, in 
an ON-OFF pattern of 1 minute on target followed by 1 minute 
on sky, with a throw of 32". Flux calibration was carried out us- 
ing Mars, and local pointing and secondary flux calibration was 
carried out using nearby bright quasars. The data were reduced 
using the facility reduction software, MOPSIC^] We list the 1 .2 
mm continuum fluxes in Table IB. 21 



Appendix C: Data collected from the literature 



In the Tables C.l C.2 C.3 and C.4 we list the X-ray luminosi 



ties, millimeter continuum fluxes and CO line fluxes that we col- 
lected from the literature, as well as their references. 



http://www.iram.es/IRAMES/mainWiki/CookbookMopsic 
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Table C.3. I2 CO 2-1 line fluxes collected from the literature; 
references listed: (1) Oberg et al. d20T0t , (2) Is ella et al. d20T0l l, 
(3) Panic et al. (2009), (4) Hughes et al. d2008T l. 



Star 


l2 CO 2-1 


Error 


Instrument 


ref. 




(Jy km/s) 


(Jy km/s) 






HD 31648 


22.0 


0.2 


SMA 


(1) 


HD 36112 


12.9 


2.58 


SMA 


(2) 


CQTau 


3.10 


0.18 


SMA 


(1) 


HD 135344B 


10.39 


0.21 


SMA 


(1) 


HD 139614 


7.00 


2.11 


RxA3/JCMT 


(3) 


HD 142527 


20.76 


0.23 


SMA 


(1) 


HD 142666 


3.12 


0.5 


SMA 


(22) 


HD 163296 


6.90 




RxA3/JCMT 


(3) 


HD 169142 


21.30 




RxA3/JCMT 


(3) 


49 Ceti 


2.0 


0.3 


SMA 


(4) 



Table C.4. 12 CO 3-2 line fluxes collected from the literature; ref- 
erences listed: (1) Dent et al. (l2003T l and (2) Panic et al. J2009l . 
(a): JCMT flux dominated by emission from the dark cloud. 



Star 


12 CO 3-2 


Error 


Instrument 


ref. 




(Jy km/s) 


(Jy km/s) 






AB Aur 


143.20° 


3.29 


RxB3/JCMT 


(1) 


HD 31648 


53.00 


1.16 


RxB3/JCMT 


(1) 


HD 35187 


< 5.00 




RxB3/JCMT 


(1) 


HD 36112 


15.70 


1.93 


RxB3/JCMT 


(1) 


CQTau 


6.00 


0.97 


RxB3/JCMT 


(1) 


HD 100546 


178.4 


19.20 


APEX 


(2) 


HD 135344B 


18.80 


0.77 


RxB3/JCMT 


(1) 


HD 139614 


9.10 


2.12 


RxB3/JCMT 


(1) 


HD 141569A 


14.70 


1.16 


RxB3/JCMT 


(1) 


HD 142666 


13.90 


2.70 


RxB3/JCMT 


(1) 


HD 150193 


1.40 




RxB3/JCMT 


(1) 


HD 163296 


8.30 


1.93 


RxB3/JCMT 


(1) 


HD 169142 


32.9 


2.51 


RxB3/JCMT 


(1) 


49 Cet 


6.58 


1.35 


RxB3/JCMT 


(1) 
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